Net current generation in a ID quantum ring at zero magnetic field 



O 
O 
(N 

X) 



I 

O 

o 



> 
oo 

o 
o 



o 



X 



Sigridur Sif Gylfadottir,^ Marian Ni^a,^ Vidar Gudmundsson/ and Andrei Manolescu^ 

^Science Institute, University of Iceland, Dunhaga 3, IS- 107 Reykjavik, Iceland 
^Institute of Physics and Technology of Materials, P.O. Box MG7, Bucharest- Magurele, Romania 

We study a non-adiabatic excitation of an electron system in a ID quantum ring radiated by a 
short THz pulse. The response of two models, a continuous and discrete, is explored. By introducing 
a spatial asymmetry in the external perturbation a net current can be generated in the ring at a zero 
magnetic field. Effect of impurities and ratchets are investigated in combination with symmetric 
and asymmetric external excitation. 



I. INTRODUCTION 

Non-adiabatic generation of current in coherent quan- 
tum rings on the nanometer scale has attracted atten- 
tion recentlypiiSi^ A one-dimensional (ID) quantum ring 
in its ground state exhibits a circulating current (persis- 
tent current) if the ring is placed in an external magnetic 
field perpendicular to the plane of the ring. The magnetic 
field breaks the left/right symmetry and the occupation 
of states carrying equal but opposite current is distorted. 
By subjecting the ring to a strong perturbation the per- 
sistent current can be changed non-adiabatically as the 
occupation of the single-electron states changesi*^. 

Here we consider two models of a one-dimensional ring, 
continuous and discrete, and investigate the response to 
a short-lived THz perturbation. The electrons confined 
to the ring are spinless and Coulomb interactions are ne- 
glected. The main objective is to see whether a current 
can be induced in the ring at a zero magnetic field by in- 
troducing a left/right asymmetry in the radiation pulse. 
The comparison between the two models shows identical 
results at least for the system parameters considered in 
the present work. We will also look at how impurities 
and ratchet-like potentials affect the resulting current. 



II. THE MODEL 

A. Continuous Model 

We consider a ID quantum ring of radius tq contain- 
ing a few noninteracting, spinless electrons. The ring is 
described by the following unperturbed Hamiltonian in 
polar coordinates: 



with I = 0, ±1,±2, .... The eigenvalue equations for the 
two operators are 



L,MO) = -iHi{e). 
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Except for ^ = all the states are doubly degenerate and 
carry a net persistent current that is normally quantified 
by the expectation value of the orbital magnetization op- 
erator 

Mo = ^vx] = -^L,. (5) 
Zc h 



B. Discrete Model 

The discrete Hamiltonian of the ID ring^ is 

Ho= 2VY,\n)H 

n 

-Vj2[\n){n+l\ + \n){n-l\], (6) 

n 

and the angular momentum operator 

L. = ^^J2^\n){n + l\-\n){n~l\]. (7) 

n 

The vectors \n) with n = 1 ■ ■ ■ N define the N points 
of the ring of radius tq separated by a linear distance 
a ^ 2TTro/N. The energy unit is V = ji^ma^). The 
common set of eigenvectors of the two commuting oper- 
ators, Hq and Lz, is 
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where m* is the effective mass of electrons. The Hamil- 
tonian is symmetric under any rotation transformation. 
Consequently Hq commutes with the angular momentum 
operator = —ihds. The common set of the eigenfunc- 
tions for the two commuting operators i?o and are 



-iW 



\J2-KrQ 



(2) 



with / = 0, ±1, ±2, ±(7V/2 - 1), 7V/2 (for N even) and 
= ^jp- is the angle corresponding to the point n of the 
ring (in polar coordinates) . The eigenvalue equations for 
the two operators are 

Ho\M = Ei\cpi) = {2V - 2Vcos{2ttI/N)) (9) 
UN 

i.lV'/) =-^sin(27r//A^)|Vi). (10) 
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In the limit l/N ^ the eigenvalue set from the pre- 
vious equations recover the values from the continuous 
model (eqs. (pj and 1^1 )■ If the ring contains impurities 
or ratchets, its Hamiltonian will be completed with the 
diagonal energies e„ 

Ho^ Ho + Y,^n\n){n\ (11) 

n 

at chosen points on the ring. We define impurities 
by adding a diagonal energy e„ = e^^p at three se- 
lected points n = 1,20,40. In case of ratchets^, four 
equidistant saw-tooth potentials are placed at the points 
n — 5, 15, 25, 35 of the ring. The saw-tooth potential 
is defined by a series of three consecutive energies: 

e„-i = O.lert, Cn = 0.2ert, (n+i = O.Sert. (12) 

The net current in the discrete model is also calculated 
from the expectation value of the orbital magnetization 
(eq. ®). 

C. Time dependent perturbation 

At t = the quantum ring is radiated by a short THz 
pulse of duration tf ~ lO/F and frequencies uji, uj2- 

Hr{t) = V{9) W{t) (13) 
W{t) = e"" sin(wii) sin(w2t). (14) 

The spatial component of the external pulse V{0) is given 
by a combination of a dipole and a rotated quadrupole 
with amplitude A: 

V{9)=A{cose + cos2{e + (t)o)). (15) 

The defiection angle (j)o between the dipole and 
quadrupole fields makes the external perturbation asym- 
metric for any values (fio ^ 0, tt. The perturbation is 
shown in Fig. To follow the time-evolution of the sys- 
tem we use the density matrix to specify the state of the 
system. The ground-state density operator p(t — 0) is 
constructed in terms of a given set of the eigenvectors 
I a) of the initial time independent Hamiltonian 

p(t = 0) = ^/(e„-Ai)|«)(a|, (16) 

a 

with the Fermi distribution /. The equation of motion 
for the density operator 

indtp{t) = [Ho + Hr{t),p{t)] (17) 

is inconvenient for numerical evaluation so we resort 
instead to the time-evolution operator T, defined by 
p{t) = T{t)p{0)T''{t), leading to a simpler equation of 
motion 

ihdtT{t) = H{t)T{t) (18) 
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FIG. 1: (Color online) Spatial distribution of the external 
potential pulse as a function of time for (a) </>o = and (b) 
(j)o — 45°. The parameter values are listed in the caption of 
Fig. El Bright /yellow regions indicate maxima in the poten- 
tial. 

with the initial condition T(0) — 1. We solve the equa- 
tion in the truncated basis of eigenvectors V'; (eq. Q) 
for the continuous Hamiltonian (eq. (Ql) or in the finite 
basis of eigenvectors of the discrete Hamiltonian (eqs. © 
or mill l using the Crank-Nicholson schemed We calcu- 
late the time evolution of the system and extend espe- 
cially the calculation of the induced current in the ring 
to times after the dependent perturbation vanishes. For 
quantifying the currents induced by Hr(t) we calculate 
the expectation value of the orbital magnetization Mo 
(eq. m in terms of the density matrix 

Mo{t)^Tr[Mo- p{t)]. (19) 

When the perturbation vanishes the time evolution of the 
system will be described by the time-evolution operator 
generated by the non-perturbed (and time independent) 
Hamiltonian Hq. The energy of the system will be con- 
stant 

E{t > tf) = Tr[Fo • p{t)] Eipuitf), (20) 

and also the magnetization for a ID pure ring 

Mo{t > tf) = Tr [Mo • p{t)] = -PbEi Ipuitf)- (21) 

We note that the non-zero values of the induced mag- 
netization in Eq. lf2T|l can be obtained only by different 
occupation of opposite current carrying states tpi and ip-i 
{pu 7^ P~i-i at i ~ t^). On the other hand, the magneti- 
zation for a ring with a finite width will oscillate in time 
in the presence of a perpendicular magnetic field due to 
coupling of radial and angular modes of density oscilla- 
tions, plasma oscillationsii 

III. RESULTS 

We use the effective mass of an electron in a GaAs, 
TO* = 0.067TOe. For modeling the ring with = 3 elec- 
trons we set the radius to rp = 14nm and in the discrete 
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FIG. 2: (Color online) The time evolution of the orbital mag- 
netization A4 o{t) calculated with the discrete model (upper) 
and continuous model (lower) for the ID ring described in 
the text. Three values of the deflection angle are presented, 
(j>o = 0°,20°,45°. The parameters of the external pertur- 
bation are huj-i = 2.83meV, ^1^2 ~ S.llmeV, KT — Ahtoi = 
11.32meV and A = 67.68meV. The frequencies are tuned to 
be comparable with the Bohr frequencies ?!.aJo,i = 2.89meV 
and Tiuii^2 = 8.6meV. The inset presents the time evolution 
of the absorbed energy (in meV) of the discrete system for 
cj>o = 0° , 45° and the time dependent part of the external 
perturbation W(t). 



model use A'' = 40 points. Initially the system is in its 
ground state where right and left rotating states, | / ) and 
I — Z), are equally occupied and the current consequently 
zero (A^o(O) = 0). In Fig. [21 we see that if the external 
pulse applied to the ring, like the ring itself, is left / right 
symmetric (0o = 0, tt) a current can not be induced since 
the occupation of states remains symmetric. However, by 
introducing an asymmetry to the perturbation a constant 
current along the ring appears in both the continuous and 
discrete models. 

For a weak excitation a perturbation analysis to 2nd 
order (in no external magnetic field) reveals that, for a 
single electron initially in the ground state | ) , the dif- 
ference between the transition probabilities to the states 
1 1 ) and I — 1 ) is proportional to sin20o 



VQs{t)-Vo,-i{t) = F{t) sin2(/.o 



(22) 



where F{t) is a double integral depending on the param- 
eters of the external perturbation. This implies that for 
a deflection angle of 0o = 45° the resulting current is at 
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FIG. 3: The induced values of the orbital magnetization Mo 
vs. the deflection angle 4>o for three different values of V (con- 
tinuous model). The external perturbation parameters are 
the same as in Fig. 2. The scale of the two lower curves 
has been magnified as indicated. The inset depicts the time 
dependent part of the external perturbation for the three sit- 
uations. 
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FIG. 4: (Color online) (a) The orbital magnetization Mo{t) 
as a function of time calculated with the discrete model of the 
ID ring. The external perturbation is asymmetric {4>o ~ 45°) 
with the same parameters as in Fig. 2 (except for A=135meV). 
The ring contains four ratchets with ert = V (explained in 
text). The almost constant blue line is the time evolution of 
Mo{t) for a pure ring, (b) The ring with ratchets plus impu- 
rities present (e™p = 0.2V). V = h^/2m*a^ = 117.5meV. 



maximum. Fig.|2|shows the induced magnetization (con- 
tinuous model) after the perturbation vanishes {t > tf) 
as a function of the deflection angle (po for three values of 
r. Clearly for strong excitation (lower T in Fig. the 
results do not depend on 0o as eq. lf22|l shows, but have a 
more complex dependence on the exact form of the per- 
turbation. As r decreases the maximum of the current 
shifts away from (/)o = 45° . The identical behavior is seen 
in the discrete model. 

The constant magnitude of the current generated in 
the system by a left/right asymmetric excitation can be 
destroyed by the inclusion of ratchets and impurity po- 
tential as can be seen in Fig. 01 In Fig. Owe see that the 
frequency of current oscillations depends heavily upon 
the strength of the ratchet potential. As the height of the 
ratchets increases the frequency grows. Furthermore the 
reduction of the current occurs even though the ratchet 
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FIG. 5: (Color online) The orbital magnetization A4o{t) as a 
function of time calculated with the continuous model of the 
ID ring. The external perturbation is asymmetric {(po = 45°) 
with the same parameters as in Fig. 2. The ring contains 
four ratchets (explained in text). The strength of ratchet 
potential trt = Eq in (a) and ert = 3.33i5o in (b). The dashed 
red curve is the time evolution of A4o{t) for a pure ring. Eo = 
h'^/2m''rfi = 2.90meV. 
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FIG. 6: (Color online) Oscillations of the orbital magnetiza- 
tion A4o{t) calculated with the discrete model. The ring is 
radiated with a symmetric pulse {(j)o = 0°). The perturbation 
parameters are the same as in Fig. 2 (except for yl=135meV). 
(a) Four sticks symmetrically placed on the ring. Each stick 
is defined by adding a diagonal = 0.2V at points 

n = 5, 15, 25, 35 of the ring. The constant line is the time 
evolution of A4o{t) for a pure ring, (b) Four ratchets on 
the ring (e^t ~ V). (c) Same as (b) plus impurities present 
(e«mp = 0.2V). V = h^l2m'a^. 



potential is very weak. It is interesting that the period 
of oscillation can be much larger than that of the applied 
perturbation. The reader may now wonder if it would 
be possible to excite a current in the system by apply- 
ing a left/right symmetric external electric field, but at 
the same time having an asymmetric static potential on 
the ring. We have tried several asymmetric potentials on 
the ring and different strength of the external excitation, 
but have not been successful in generating a steady DC 
current. Some results are seen in Fig. El Clearly, as we 
have drawn attention to above, the frequency of the cur- 
rent oscillations can be made very low compared to the 
dominant frequencies included in the excitation pulse. 

IV. CONCLUSIONS 

For a quantum ring in no external magnetic field a 
left/right asymmetric excitation is essential for the ap- 
pearance of a circulating current. In an external mag- 
netic field the energy of states whose orbital angular mo- 
mentum is aligned with the field is lowered, thus break- 
ing the symmetry and generating a DC current. Here we 
have seen that a similar effect can be accomplished by an 
asymmetric excitation pulse at zero magnetic field, but a 
DC current can not be generated by a left/right symmet- 
ric excitation even though the ring itself is made asym- 
metric by a static potential or by change to the geome- 
try, breaking the perfect circular symmetry. The intro- 
duction of a ratchet potential in the ring that is excited 
by an asymmetric external field changes the DC current 
created into an AC current with arbitrary low frequency 
depending on the intensity of the ratchet. This is a clear 
nonlinear behavior that could not be observed by tradi- 
tional linear response methods. 

The electron-electron interaction has been neglected 
here as the effects studied do not depend on it qualita- 
tively, but certainly quantitatively they doi^ 
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